Vortex state in double transition superconductors 
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Novel vortex phase and nature of double transition field are investigated by two-component 
Ginzburg-Landau theory in a situation where fourfold-twofold symmetric superconducting double 
transition occurs. The deformation from 60° triangular vortex lattice and a possibility of the vortex 
sheet structure are discussed. In the presence of the gradient coupling, the transition changes to 
a crossover at finite field. These characters are important to identify the multiple superconducting 
phase in PrOs4Sbi2. 
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In the newly discovered superconductor PrOs4Sbi2, 
which is a heavy-fermion compound with filled skutteru- 
dite structure, a novel pairing mechanism is considered 
in connection to quadrupole fluctuations [ij. Several ex- 
perimental results suggest unconventional pairing sym- 
metry. In the nuclear spin-lattice relaxation experiment, 
the coherence peak is absent below the transition tem- 
perature T c , and the low temperature behavior indicates 
that PrOs4Sbi2 has full gap or point nodes, excluding line 
nodes, at a zero field [3- The specific heat jump at T c 
shows double transition |2| . The thermal transport mea- 
surement in the magnetic field rotated in the a6-plane of 
the crystal axes reveals the phase diagram of the dou- 
ble transition, and shows that a fourfold-symmetric pair- 
ing function around the c-axis in the high-field H-phase 
is changed to a twofold-symmetric one in the low-field 
L-phase at the second transition field H* , when mag- 
netic field or temperature is decreased 0. Thus, the 
H-T phase diagram is divided into two regions, H- and 
L-phases. The muon spin relaxation (/xSR) experiment 
reports the spontaneous moment in the superconducting 
phase 0| . These results indicate an unconventional pair- 
ing of the superconductivity in PrOs4Sbi2. 

To explain the double transition, we need multiple 
components for the pairing functions. The fourfold- 
twofold transition in PrOs4Sbi2 can be explained as 
follows: After the fourfold-symmetric pairing compo- 
nent appears at the first transition, the second compo- 
nent appears at the second transition, and the combi- 
nation of these components gives rise to the twofold- 
symmetric gap structure. For PrOs4Sbi2, the scenario 
of "anisotropic-s"+id-wave pairing with point nodes was 
proposed, and the double transition is analyzed by the 
Ginzburg-Landau (GL) theory at a zero field 0. In 
the vortex states at finite fields, however, this fourfold- 
twofold transition is not a true transition, when the gra- 
dient coupling terms are present in the GL theory Q. 
That is, the twofold component at the L-phase is in- 
duced by the gradient coupling, and produces twofold 
symmetric behavior in addition to fourfold symmetric 
one, even in the H-phase above H* . To avoid the mix- 
ing of the twofold component above H* , the plausible 



pairing symmetry should be a triplet pairing where d- 
vector of the two components are orthogonal each other 
(i.e., dt - d 2 = 0), so that the gradient coupling terms 
vanish y|. 

On the other hand, in multi-component superconduc- 
tors, it is interesting to examine the possibility of exotic 
vortex states, such as a coreless vortex 0, 0, 0, 0, 0, 
In this letter, we investigate the double transition 
by the two-component GL theory in the magnetic field, 
and calculate the vortex structure by the time-evolution 
of the time-dependent GL (TDGL) equation. 

First, we derive two-component GL equation appropri- 
ate to study the fourfold- twofold transition in PrOs4Sbi2 . 
The pair potential in 2 x 2 matrix form is decomposed 
to two components as A(r,k) = r]i(r)0i(k) + 7/ 2 (r)(/>2(k) 
with the order parameter ry m (r), where r is the center 
of mass coordinate of the Cooper pair and m — 1,2. 
The relative momentum k of the pair is mapped on 
the Fermi surface. The pairing function is given by 
4> m (k) — ia y 4> m (\i) for the singlet pairing, and <^ m (k) = 
i 2 J=a . z dmjQtjaj&y for the triplet pairing with Pauli 
a y , a z . We assume that the superconduct- 
mp Dy </>i(k) (</> 2 (k)) has fourfold (twofold) symme- 
try [a, |tJ , and that the transition temperature estimated 
from the pairing interaction is lower for the second com- 
ponent, i.e., T c — T c i > T C 2. Since the pairing symmetry 
for PrOs4Sbi2 is not established yet, the pairing function 
forms </> m (k) are not specified in this study. 

Within the GL approximation, the free energy in the 
superconducting state is generally given by F s — F n + 
J /(r)dr with 

/(r) = -a (T c - T)\ m \ 2 - a (T c2 - T)\ m \ 2 



matrices a x 
ing gap by 
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in the clean limit 0], where q = (7i/i)V + (27r/</>o)A, 
F n is the free energy in the normal state, \A%\ = 
7C(3)/(167r 2 T c 2 ) with Riemann's ^-function, <fio is a flux 
quantum, v is a Fermi velocity, (• • •) indicates the Fermi 
surface average of k, and A is a vector potential. Since 
the magnetic field is applied along the z-axis, q z = 0. 
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In the dimensionless form, Eq. QJ is written as 

/ (, TV „ 



+012^^772 + C 12x r) 2 q x rn 

+Ci2 V vi<ilm + CwyvZqlm 

+l{\vi\ A + C 2 \V2\ 4 + 4C 3 \ m \ 2 \ m \ 2 
+Ctr,?r,l + C 4 7/ 2 * 2 t/ 2 }, 



(2) 



using /o = aoT^l rj 2 = a T c /(2|A 2 |C 1 ) = 1 and £ 2 = 
M^ICiiz/^oTc) = 1. The coefficients of the gradient 
terms in Eq. (J2J are given by 
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where 1 



/>^0n for the singlet pairing, 
12 



mrn — u m 

The slope 



d„ for the triplet pairing and \(f) m \ = v-m'+'m 
ratio of H* and H C 2 in the H-T phase diagram is roughly 

given by X = [(^TF^)(WF«. 2 )/(WP^)(WF^)]; /2 - 
To reproduce the experimentally obtained H* behavior, 
it is appropriate to use X ~ H* / H C 2 ~ 0.5. In Eq. 
©, c is an anisotropic parameter related to the second 
order parameter. When |</> 2 | 2 has two- fold symmetry un- 
der the rotation around the z-axis, c is not zero since 
)■ C\2x and C\2 V reflect the strength 
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of the gradient coupling between 771 and 773- 

The coefficients of the quadratic terms in Eq. 
are given by & = (|0i| 4 ), C 2 = (|0 2 | 4 }/C 1 , C 3 = 
(I0i| 2 |0 2 | 2 )/Ci, C 4 = (0f 0l)/Ci for the singlet pairing, 
andCi = (2|d 1 | 4 -|d r d 1 | 2 },C 2 = (2|d2| 4 -|d 2 -d 2 | 2 }M, 
C 3 = (|d 1 | 2 |d 2 | 2 + \d{ ■ d 2 | 2 - |d x • d 2 | 2 )/d, C 4 = 
(2(d?| • di) 2 - (d^ • d|)(di • di))/Ci for the triplet pair- 
ing. However, we can not identify the definitive values 
of coefficients in Eq. @ for PrOs 4 Sbi2, because the de- 
tailed information of the pairing function and Fermi sur- 
face structure have not been established yet. It is noted 
that the Fermi surface anisotropy also largely affects on 
the coefficients. Therefore, we treat these coefficients as 
arbitrary parameters, and report some typical results ob- 
tained in this framework. 

Before considering vortex states, we study a uniform 
state at a zero field. From the free energy minimum 
condition, the relative phase of 7/1 and 772 becomes (a ± 
7r)/2 with a given by C4 = |C4|e la . For the H-phase, 
7/1 = (I — T/Tc) 1 / 2 and 7/2 = 0. The second component 
772 appears at the lower transition temperature T* given 
by 



T* = T c2 /T c -(2C 3 -\C 4 \) 
T c 1 - (2C 3 - |C 4 |) ' 



(4) 



which is derived by linearizing the equation df /drfe = 0. 
In the presence of 771 , T* is suppressed compared with 
T c2 . To assure that T* > 0, we have to satisfy 2C 3 — 
\Ci\< T c2 /T c . 

The vortex structure is calculated by the time- 
evolution following the TDGL equation coupled with 
Maxwell equation 
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—A = j s - k 2 V x B, 



B = V x A. 



(5) 
(6) 



The supercurrent j s 
given by 



(js,xJa,v) K (df/9A x ,df/dA y ) 



j SfX = Re[r)*(q x rn) + C^x^Hqxil-i) 

+c 12 xvl(qxV2) + c^xvUixm)}, (?) 

3b, v = Re [r)$(q y T]i) + C 2 2yV2(<lyV2) 

+c 12y iil(qym) + CwyvZiivVi)]- ( 8 ) 

We use the same scale units as in Refs. and 
for length, field and time, except for the order param- 
eters. We here scale rj m by 770 instead of r] (T) = 
7/o(l — T/Tc) 1 / 2 . In our calculations, we typically use 
the GL parameter k, = 4, which belongs to a high-K case, 
i.e., the order parameter structure is not significantly af- 
fected by the internal field distribution. 

Calculations are performed in a two-dimensional 
square area. Outside the open boundary, we set 7/1 = 
7/2 = and B(r) = H with an applied field H. For the 
initial state of 771, 772 and A inside, we use a uniform state 
at a zero-field. Through the time evolution, vortices pen- 
etrate from the boundary. After enough time later, the 
time-evolution converges and the vortex lattice state is 
obtained. We analyze this final state. 
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FIG. 1: (a) Maximum of I772 (r) | (solid lines) as a function of 
H/H c2 {0). T/T c = 0.1(o), 0.4(«) and 0.7(A). 7/2 appears at 
H < H*(T). For T/T c — 0.1, H-dependence of |r7i| is also 
presented (a dashed line), (b) H-T phase diagram in this GL 
theory. The transition field H* and H C 2 are presented. 

We discuss the case when the gradient coupling is ab- 
sent (Ci2x = Ci 2y = 0). Here, we set T c2 /T c = 0.95, 
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C 2 = 1,C 3 = -G\ = 0.2, c = 0.3 and X = 0.5. In Fig. 
I2a), we show the maximum order parameter amplitude 
at T/T c = 0.1, 0.4 and 0.7 as a function of the external 
field H/H c2 (0) with H c2 (0) = <?V(27r£§). At a critical 
field H*(T), the second component 772 appears as a sec- 
ond order transition. The T-dependence of H* is shown 
in Fig. [Jb). This H-T phase diagram qualitatively re- 
produces the double transition of PrOs4Sbi2 |4|. 
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FIG. 2: (a) Density plot of an internal field distribution B(r) 
at H/Hc2(0)=0.2 and T/T c = 0.1. White region corresponds 
to the vortex core with large B(r). Density plot of |??i(r)| (b) 
and \r/2 (r) | (c) around a vortex. Diffraction pattern |S q | 2 at 
H/H c2 {0)=0.2 [L-phase] (d) and 0.5 [H-phase] (e). 

At H/H c2 {0) = 0.2 and T/T c =0.1, the vortex dis- 
tribution is presented in Fig. Efa), where the internal 
field distribution is shown as a density plot. There are 
some domains of the different orientation of the trian- 
gular lattice, since the vortex lattice configuration is af- 
fected by the boundary. And the triangular lattice is 
deformed from 60° triangle by the effect of finite c in Eq. 
(|3*|) , coming from the twofold symmetry of 4>2 ■ The vor- 
tex core shapes of |r?i(r)| and 1 772 (r ) | are, respectively, 
shown in Figs. Efb) and[2c). Vortex core shape of 
\f] 2 (r) I is stretched out toward the y-direction, while that 
of 1 771 (r) I remains to be circular. The deformation of 
the vortex lattice is clear, when we see the form factor 
of the internal field distribution, as in the neutron scat- 
tering experiments. Figure |2td) shows the "diffraction 
pattern" |i? q | 2 in the L-phase, with the Fourier compo- 
nent Bq = ^ IT e 1<1 ' r B(r), where we exclude the region 
near the boundary in the integration over r. There ap- 
pear six peaks of the triangular lattice, and the rotated 
ones for different orientations. These peaks are on an el- 



lipse with long axis 2q x and short axis 2q y with the ratio 
Qx/% ~ 1-2. For a reference, we also show the diffraction 
pattern of the 60° triangular lattice in the H-phase in Fig. 
|2te), where peaks appear on a circle. The observation of 
this vortex lattice deformation may be another means to 
detect the fourfold-twofold transition in the H-T phase 
diagram. 
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FIG. 3: Vortex structure at H/H c2 (0) = 0.3 and T/T c = 
0.1 [L-phase]. Inside region is focused, (a) Density plot of 
1 771 (r) I and 1 772 (r) | - In the black regions, |?ji| and |?72| share 
a vortex core. Green (red) regions show that only \rji\ (\r]2\) 
has a vortex core, (b) Internal field distribution B(r). (c) 
Relative phase arg{7/i(r)/?72(r)}. Right panel schematically 
shows relative phase near the vortex sheet appearing at the 
domain wall. 

At higher field in the L-phase, we also find the vortex 
sheet structure Q in addition to the regular vortices. 
We show the spatial distribution of 1 771 (r) j and (772 ( r )| 
at H/H c2 (0) = 0.3 in Fig. da). The black circle re- 
gion presents regular vortex, where |?7i(r)| and 1 772 (r) | 
share the same vortex core, as seen in Fig. |2Ja). The 
green (red) circle region shows the rji- (r) 2 -) vortex, where 
only |?7i(r)| ( I772 (r) | ) has vortex core and the other | r) 2 (r)\ 
(|?7i(r)|) does not. These green and red vortex cores are 
located alternatively along a loop, forming vortex sheet. 
In the internal field distribution shown in Fig. \j$b), B(r) 
has a sharp localized peak at a regular vortex. And large 
B(r) regions of the vortex core on the vortex sheet are 
connected each other along the vortex sheet. Each of r\\- 
and 772-vortices has half flux quantum. If this line struc- 
ture is found by the direct observation of the internal 
field distribution, it can be evidence of the vortex sheet 
appearing in unconventional superconductors. The rela- 
tive phase of r?i(r) and 772(1") are presented in Fig. EIc). 
Around the regular vortex, the relative phase is fixed 
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at (a + 7r)/2 or (a — ir)/2 (in our parameter, a = n). 
Across the vortex sheet, the relative phase changes from 
(red region) to it (blue region). This indicates that 
the vortex sheet appears at the domain wall between the 
region with the relative phase (a + tt)/2 and that with 
(a — 7r)/2, and that it is not easy to disconnect vortices 
on the vortex sheet. As shown in the right panel in Fig. 
|3Jc) , since windings of the relative phase are opposite 
at the 771-vortex and the 772-vortex, the relative phase 
changes from to ir through ir/2 (yellow) or 3ir/2 (pur- 
ple) alternatively between the nearest neighbor vortices 
along the vortex sheet. Since domains with relative phase 
(a ± 7r)/2 are degenerate in free energy, we can expect 
both domains to coexist in sample materials. In the pres- 
ence of the domain wall between domains, vortex sheet 
appears when applying field. In our simulation for the 
penetration process of vortices, boundary region helps 
creation of the domain wall with vortex sheet structure. 
The domain wall and vortex sheet are tight structure 
when coming inside, and survive stably. 
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FIG. 4: Maximum of 1 772 (r) | as a function of H / H C 2(Q) in the 
presence of the gradient coupling. T/T c = 0.1(o), 0.4(») and 
0.7(A). rj2 survives up to H C 2(T). 

Lastly, we consider the case when the gradient coupling 
terms do not vanish. We show the results when C\2 X — 
0.2 and C\2 y = and the other parameters are kept same. 
We note that, if we consider the s+id-state (e.g. <?!>i(k) cx 
l-(k* + k* + k 4 z ) and^k) oc A£-fc*) @, C 12y = 0. We 
see the similar vortex structure as in Figs. [2and|3]in the 
L-phase. The field dependence of the order parameter 
maximum of 1 772 (r) | is shown in Fig. 0] On raising field, 
1 772 1 decreases, but it does not vanish at H* . 1 772 1 has 
finite value up to H C 2(T). Therefore, twofold symmetric 
order parameter 77 2 ( r )02(k) is mixed in addition to the 
fourfold one 77i(r)(/>i(k) in the H-phase 7]. That is, H* is 
not a phase transition but the crossover field where 772 is 
enhanced in the vortex state. While the anomaly of the 
second order transition appears at T* in the zero field 
case, the anomaly of the phase transition is not observed 
around H* at finite fields. 

In summary, we investigate the vortex state based 
on the two-component GL theory in the situation when 
the fourfold-twofold symmetric superconducting transi- 
tion occurs. We estimate the transition field H* where 
second order parameter r\2 appears. However, H* is 



changed to a crossover field when the gradient coupling 
terms exist, because the small 772 survives up to H c2 (T). 
In the L-phase below H*, the vortex lattice deforms from 
60° triangular lattice due to the effect of the twofold sym- 
metric second order parameter. In the two-component 
superconductor, there is a possibility to observe the ex- 
otic vortex structure such as vortex sheet at the domain 
wall, where two order parameters have different vortex 
cores and these cores are alternatively located along a 
line. These characters of the vortex structure may be 
clear evidence of the fourfold-twofold symmetric double 
transition and unconventional multi-component super- 
conductivity, if they are experimentally observed. 
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